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Abstract 

We derive the Hamiltonian for spherically symmetric Lovelock gravity using 
the geometrodynamics approach pioneered by Kuchaf [1] in the context of four- 
dimensional general relativity. When written in terms of the areal radius, the 
generalized Misner-Sharp mass and their conjugate momenta, the generic Love- 
lock action and Hamiltonian take on precisely the same simple forms as in gen- 
eral relativity. This result supports the interpretation of Lovelock gravity as the 
natural higher- dimensional extension of general relativity. It also provides an 
important first step towards the study of the quantum mechanics, Hamiltonian 
thermodynamics and formation of generic Lovelock black holes. 



1 Introduction 



The elegance of its underlying geometrical structure makes Einstein's general theory of 
relativity both beautiful and compelling. More importantly, observational consequences of 
general relativity have been verified over a large range of macroscopic scales, from terrestial 
to cosmological. The inevitability of singularities in the context of cosmology and black 
holes makes it clear that modifications to the Einstein-Hilbert action are inevitable at very 
short distance scales. Quantum theory suggests the need to add higher-curvature terms 
whereas string/M-theory [2] requires the existence of extra spatial dimensions. The most 
natural such generalization of the Einstein-Hilbert action in n dimensions that is quasi- 
linear second-order is provided by so-called Lovelock gravity [3], which consists of a sum of 
[n/2] dimensionally extended Euler densities that are non-trivial only in dimensions greater 
than four, ([n/2] refers to the largest integer less than or equal to n/2.) 

The fact that Lovelock gravity is second-order implies that it is ghost-free when lin- 
earized around a flat background [3]. Moreover, string theoretic arguments [5] suggest 
that quadratic Lovelock gravity appears in the low-energy limit for strings propagating 
in curved spacetime. Since Lovelock gravity is a natural extension of general relativity in 
arbitrary dimensions, it could potentially provide insights into quantum gravity or even 
a hint concerning the origin of the four- dimensionality of our universe. Lovelock gravity 
has therefore received considerable attention recently leading to interesting developments 
at the classical level. (See [6] for review.) There are, however, few results on its quantum 
aspects. The purpose of this letter is to pave the way towards a systematic study of the 
quantization of Lovelock gravity. 

At first glance, the action for Lovelock gravity appears quite complicated [3j: 

1 f [n/2] 
n J p=0 

(p) ' 2p Pi-pp°'i- cr p ^i^i •■• 1 V P "p ' \ L - Z ') 

where SplZp* '■= ^[ p \---^pp] ■ Each is a coupling constant of dimension (length) 2 (p _1 ). 
a!(o) is proportional to the cosmological constant and the p = 1 term corresponds to general 
relativity. As stated above, the specific combination of curvature terms in each of the 
guarantees that the equations for Lovelock gravity are second order in the metric compo- 
nents and ghost-free [4J. Equally important is that the theory obeys a generalized Birkhoff 
theorem [8]: for a given set of coupling constants which do not admit degenerate vacua, 
all C 2 -class spherically symmetric vacuum solutions except for the Nariai-type solution ad- 
mit a hypersurface orthogonal Killing vector and can be parametrized by a single constant 
of motion. The spherically symmetric vacuum solutions have been fairly well studied [6], 
but little work has been done on the quantum mechanics perhaps because of the perceived 
complexity of the equations. 
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The Hamiltonian analysis for full Lovelock gravity, first considered by Teitelboim and 
Zanelli [10J, is qualitatively different from that of general relativity due to the higher-order 
nature of the theory [TOl [TTj . For the case of spherical symmetry, the geometrodynamics 
of five- dimensional Einstein- Gauss-Bonnet (i.e. p = 2 Lovelock gravity) was worked out 
by Louko et al [12], while the Hamiltonian analysis of higher-dimensional Gauss-Bonnet 
coupled to matter was recently done in [13] . 

In the following, we will show that the geometrodynamic formulation of Kuchaf [1] yields a 
Hamiltonian for spherically symmetric p-th order Lovelock gravity that is just as simple and 
elegant as that of Einstein's original theory. This result suggests that Lovelock gravity is in 
a deep sense a natural higher- dimensional extension of general relativity. It also provides 
a first important step towards a more complete understanding of the quantum properties 
of Lovelock black holes. 



2 Geometrodynamics of Lovelock Black Holes 

We start with a general diffeomorphism invariant theory of gravity that obeys the Birkhoff 
theorem in spherically symmetric spacetimes. Since an n(> 4)-dimensional spherically sym- 
metric spacetime (Ai n , g^ v ) is a warped product of an (n — 2Y-dimensional unit maximally 
symmetric space (5' n ~ 2 ,7 a &) with positive sectional curvature L 1 ! and a two-dimensional orbit 
spacetime (M 2 ,g A B), the line element may be given by 

g^dxW = g AB (y)dy A dy B + R 2 (y)dQ 2 n _ 2 , (2.1) 

where dfl 2 l _ 2 is the line element on (S n ~ 2 , jab)) an d the areal radius R is a scalar on 
(M 2 ,g AB ). In a spherically symmetric spacetime, the system may be described by the 
effective two-dimensional action; 

1(2) = J dy°L[g ABj R] = J df J dfC[ gAB ,R], (2.2) 

where y° is a timelike coordinate on (M 2 , g AB ). Here the Lagrangian L and the Lagrangian 
density £, are functionals of the metric functions and their derivatives, which are determined 
up to a total derivative. The main purpose of geometrodynamics is to find canonical 
variables (that are functionals of the metric functions) to provide a tractable form and 
transparent physical meaning for C 

If (DR) 2 := g AB (D A R)(D B R) ^ 0, where D A is the covariant derivative on (M 2 ,g AB ), the 
metric may be written in terms of the areal coordinates as 

ds 2 = -F(R, T)e 2a{R > T) dT 2 + F'^R, T)dR 2 + R 2 dtt 2 n _ 2 . (2.3) 
1 The discussion can easily be generalized to included zero or negative sectional curvature [9]. 
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The original Birkhoff theorem in general relativity asserts that the spherically symmet- 
ric vacuum solution is uniquely given by the Schwarzschild-Tangherlini solution, which is 
given in the coordinates ( 12. 3 p by F(R,T) = F(R; M) and a(R,T) = 0, where M has the 
interpretation as the mass and T is the so-called Schwarzschild time. In the presence of a 
cosmological constant, the Nariai solution is also possible where (DR) 2 = is satisfied. In 
Lovelock gravity, the most relevant part of the statement of the Birkhoff theorem is that 
the C 2 vacuum solution with (DR) 2 ^ is the Schwarzschild-Tangherlini- type solution if 
the theory does not admit degenerate vacua [7J [S] . (The theorem in the quadratic theory 
was proven in [IJJ [15].) We assume that the metric is C 2 , keeping in mind that there are 
counterexamples in the literature |16j. 

In the fully dynamical setting, we can without loss of generality consider precisely the 
same form of metric but with M = M(R, T) and o = cr(R, T) both as functions of R and 
T. In that case M(R,T) has the geometrical/physical interpretation as the generalized 
Misner-Sharp mass function. Because of (DR) 2 = F(R; M(R, T)), the mass function can 
in principle be expressed as a function of (DR) 2 , as long as (DR) 2 ^ 0, the theory which 
does not admit degenerate maximally symmetric vacua and one knows the form of the 
static spherically symmetric solution, i.e. F(R, M). In Lovelock gravity, the mass function 
is given by 

[n/2 

in — ' i -i .-) 

M 



n 2)A n -2 v-v ~ D n-i-2p ft /'nD^2^^ 



2k « p=0 

[n/2 



= (n 2)A n _ 2 ^ ~ (p)ir -i- 2 p (1 _ m m)) p ? {2A) 

n p=0 

where A n _ 2 is the area of (S n ~ 2 , 7^) and := (n — 3)\a^/(n — 1 — 2p)\ [T71 [TSJ [13 [S]- 
In n-dimensional general relativity without cosmological constant, it reduces to 

M = (n ~^ 2 An - 2 a {1) R n - 3 (1 - (DR) 2 ) , (2.5) 

and hence 

9k 2 M 

F(R, M) = 1 — fp^ -. (2.6) 

In general, the expression for F(R, M) is more complicated. In order for the theory 
to provide a unique solution for each M, the relationship between M and F must be 
monotonic, which provides a condition on the Lovelock coupling constants in ( 12. 4p . 
We will for what follows assume that this condition is satisfied, but will return to the more 
general case in [9]. 



The essence of Kuchaf ' geometrodynamics performed in vacuum four-dimensional general 
relativity is to write the Hamiltonian equations describing the dynamics of a general spher- 
ically symmetric geometry in terms of the areal radius, mass function and their conjugate 
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momenta. We now present the geometrodynamics of spherically symmetric Lovelock grav- 
ity. 



For spherically symmetric spacetimes, the action reduces to 



A _ 2 f _ [n,2] 

1 = 'YJ' J d2 vV-9M Rn * £ a (p) C (ph 



(2.7) 



where g^,) '■— det(g A B) arid the dimensionally reduced pth order Lovelock term is given 
from expressions (2.19) and (2.20) of jH] as 



C 



(p) 



,1 - (DR) 2 Y , ,D 2 R fl-(DR) 2 ^ p - 1 

(n-2p)(n-2p-l) -2p(n-2p)- ' 



i? 2 



# 2 



(n-2)! 
(n-2p)! 

+ 2p(p _ 1) {D*RY-{D a DbR){D b DaR) ( l-(DR)y - 2 | p W - ( M )>n 



i? 2 



i? 2 



i? 2 



(2.8) 



where ^72. is the Ricci scalar on (M 2 ,g A B) and D 2 R := D A D A R. The contraction was 
taken over the two-dimensional orbit space. Using the binomial expansion and collecting 
terms with the same number of derivatives yields 



(v) 



(n-2)! 
(n - 2p)\ 



( 'l-(7Ji?) 2 \ p , D 2 R {l-(DR)^ p - 1 

(n - 2p)(n - 2p - 1) v ; - 2p(n - 2p)- 1 



R 2 



+ ^ 7 " + ^(, + l)!(p-2-0!l i 



D A R{D 2 D A R - D A D 2 R){DRf 



(7J 2 i?) 2 - (D A D B R)(D B D A R) 

R 2p-2 



R 2 p- 



(2.9) 



where we have removed some terms containing ( 2 )7£ by using the two-dimensional identity: 

,( 2 ) 



(DR) 2 K = 2(7J A i?) (D 2 D A R - D A D 2 R) . 
This identity can be derived from Eq. (2.10) of 



(2.10) 



Using integration by parts, the term inside curly brackets in (12. 9 j) can be replaced by the 
following remarkably simple term: 

£ (• + ^]7JZ^ _ ^ D A (R^)((DR) 2 y^D A ((DRf) - (D A R)(D 2 R)} . (2.11) 
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Finally, the pth order Lovelock Lagrangian is, up to total divergences: 
(n-2)! 



[DRY \ 



[,,) (n-2p)! 

+ (n - 2p)(n - 2p - l)|(l - (PP) 2 ) P + 2p(PP) 2 jp~ 



2p 



(2.12) 



This is the first key result of the present paper. Its value becomes evident when considering 
the following Arnowitt-Deser-Misner (ADM) decomposition on (M 2 , qab)'- 

ds 2 = -N{t, xfdt 2 + A(t, x) 2 (dx + N r (t, x)dt) 2 + P 2 (x, t)dtt 2 n _ 2 . (2.13) 

Specifically, Eq. (I2.12p enables us to find a canonical transformation from ADM phase space 
variables (A, Pa; R, P r ) to (M, P M ; S, Ps). (Our x, S, and P$ are equivalent to Kuchaf's r, 
R and Pr, respectively.) The new canonical variable S is then chosen as S = R and hence 
also represents the areal radius. 

(DR) 2 contains no derivatives of A, while the Ricci scalar term in the action takes the form: 
, (2) r 1 . 

^^R n ' 2p n = —2N < (R n ~ 2p y - N r (R n - 2p y Ua - n',a - jv p a') 

^NUR^y'A- 1 + (ir*-*)'^- 1 )' }+$(...) + d x {...), (2.14) 

where a dot and a prime denote the partial derivative with respect to t and x, respectively. 
Thus it is evident that the generic Lovelock term is linear in A, a fact that plays a crucial 
role in the following analysis. We therefore write the effective two-dimensional Lagrangian 
density £ in (12.21) as 

C = B (R, y, A) A + B 1 (R, y, A), (2.15) 
where we have replaced R in the above by 

y:=i(p-Ar r P') (2.16) 

because general covariance of the action guarantees that R always appears in this combi- 
nation [I]. Equation ( 12.151) implies that the conjugate to A is 

P A = B (R,y,A), (2.17) 

which implicitly determines y = y(A, P\, R) as a function of (A,Pa,P). It is important 
for the following that y is independent of Pr, but the explicit form of B>o(R,y, A) is not 
crucial. 
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Following pQ, we define a new set of variables (M, P M , S, P$) from the ADM phase space 
variables (A, Pa, R, Pr) by 

P M ■= -e°T' = -^, (2.18) 

r 

S := R, (2.19) 
P s := P S (A,P A ,R) = ^(P R R' -AP A -P M M'), (2.20) 
while M is given as a function of (A, Pa, R) by solving Eq. (12.41) and 

F(M,R) = ^-y 2 . (2.21) 
We recall y = y(A, P A , R) as implied by ( 127171) . 

An explicit calculation reveals that this transformation preserves the Liouville form up to 
a total variation and is therefore canonical: 

P A 5A + P R 5R = P M 8M + P S 5S + <?(...). (2.22) 

As argued by Kuchaf [1] the specific form of P$ is determined by the transformation 
properties of the variables under spatial diffeomorphisms. Equations ( I2.18l) - (12.21l) are 
completely generic in form. The theory specific information is contained in the functions 
F(M, R) and y(A,P A ,R). 



To proceed, we note the key identity: 
D A (R n - 2 P)D A ((DR) 2 ) 



{1 - (1 - (PP) 2 )"- 1 } 



{DR) 2 

= (^-2p)^^{l-(l-Pr 1 }{-AyP+^ + iV r A^P'}. (2.23) 

By substituting (g53J) and (12TT4"j) into (I2~T2|) one obtains: 
n-2)\ 



£ = 111 - pN^A-'R 1 '^ 

(p) (n-2p-l)V 



2y(N' r A + N r A') 



9 N 

_ J (^-^/'A-i + ( j R«-2p)'( A -i)' 



(n - 2p)R n ~ 2 P 

+ P^jl - (1 - P) p ~ 1 }(AiV r2/ + A-'NR^F' - \2yA + Ay^J +(1 - F^Ay^ 

(2.24) 

Eqs. (J2HJ), ( j23HD and d2~2T|) imply that: 

n p=0 I J 

(2.25) 
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The result for PmM' is analoguous with dots replaced by primes. 



Using the above expressions for PmM and PmM' as well as (I2.24p it is straightforward, 
but algebraically lengthy, to confirm the following via direct substition: 

C- P M M -—M' = — (\ r ■ X—) M' 
M B! F \ R'J 

= ~§PmM>. (2.26) 



From Eq. (12. 26p it follows immediately that the Hamiltonian in the equivalent two-dimensional 
theory takes the form: 

U G ■= P M M + P S S-C 

= N M M' + N S P S , (2.27) 

where we have used Eq. (I2.26P to get the second line and defined new Lagrange multipliers 
N M and N s . The Lagrangian density for the canonical coordinates (M, S, N , N s ) can 
now be written as 

C = P M M + P S S - N M M' - N S P S . (2.28) 

Eq. (I2.28P is precisely the same form as Eq. (122) of [I]. In comparison to the rather 
complicated starting point in Eq. (12.81) . this equivalent Lagrangian density is extremely 
simple and the physical meaning of the canonical variables are very clear. This is our 
main result. Remarkably, the coupling constants a>( p ) do not appear explicitly in any of the 
equations after (12.181) . They are in fact hidden in the definition of the mass function. This 
makes it possible to treat any class of Lovelock gravity in exactly the same way. 



The constraints, M' = and P$ = 0, are obtained by varying the Lagrange multipliers N M 
and N s , respectively. On the constraint surface M = m(t), as expected, and Ps = 0. The 
reduced phase space is therefore two-dimensional consisting of p m := dxPM(x,t) and 
m. With suitable boundary conditions [H], one can repeat the analysis of pQ for spacelike 
slicings that intersect both left and right branches of the outer horizons of eternal black 
holes to obtain the reduced action: 



'(2) 



dt 



p m m — (N + — N^)m 



(2.29) 



where N± := ^lim x ^.± oa N M . The reduced equations of motion in vacuum then imply 
that m = m = constant, and p m = —(N + — AL). 
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3 Conclusions 



In summary, we have presented the results of a geometrodynamics/Hamiltonian analysis 
of spherically symmetric Lovelock gravity. The equivalent two-dimensional Lagrangian 
density f!2.28|) and Hamiltonian (I2.27P written in terms of the areal radius, the generalized 
Misner-Sharp mass, and their conjugate momenta have exactly the same simple forms as in 
general relativity. Our analysis paves the way for further understanding of generic Lovelock 
black holes. Details of the proof will be presented in the subsequent paper [9]. 

Given (12 .29 p . the reduced quantization of Lovelock black holes can be performed exactly 
as in [1] . Alternatively one can take the approach of Louko and Makela [19] to derive the 
quantization of the throat of the Lovelock black hole to obtain the mass spectrum as a 
function of the number of spacetime dimensions and Lovelock couplings. This is currently 
in progress. As well, the geometrodynamics described above is the starting point for the 
Hamiltonian thermodynamics of Lovelock black holes, following [12]. Finally, we remark 
that, as we will show in [9], the Hamiltonian (I2.27P allows one to write down the reduced 
Hamiltonian equations for black-hole formation via scalar-field collapse, either in gener- 
alized Painleve-Gullstrand coordinates, as done for Einstein-Gauss-Bonnet gravity in [13] 
or using the local Hamiltonian of [201 EI|- This in turn provides the starting point for a 
detailed numerical study of black-hole formation in Lovelock gravity. 
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